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Measurements of the proton structure function F2 in the Q 2 range 0.6 — 17 GeV 2 from ZEUS 1995 shifted vertex data and 
Q 2 ~ 1.5 - 20000 GeV 2 from 1996 and 1997 ZEUS data are presented. From the former and other ZEUS F 2 data the slopes 
dF^ldhxQ 2 at fixed x and dF? / d\n{\ / x) at fixed Q 2 are derived. F2 data at Q 2 below 0.9 GeV 2 are described successfully by a 
combination of generalised vector meson dominance and Regge theory. Using a NLO QCD fit the gluon density in the proton is 
extracted in the range 3 X 10~ 5 < x < 0.7 from ZEUS 1994 and 1995 data. For Q 2 ~ 1 GeV 2 it is found that the qq sea distribution 
is still rising at small x whereas the gluon distribution id strongly suppressed. It is shown that these observations may be understood 
from the behaviour of the F2 and flt/'a/dln data themselves. 



1 Introduction 

Measurements of the low and medium Q 2 neutral cur- 
rent (NC) deep inelastic scattering (DIS) cross sections 
at HERA have revealed the rapid rise of the proton struc- 
ture function F2 as Bjorken-x decreases below 10~ 2 . At 
low Q 2 down to 0.1 GeV 2 ZEUS data allows study of the 
'transition region' as Q 2 — * in which perturbative QCD 
(pQCD) must break down. At high Q 2 , NC DIS mea- 
surements are sensitive to details of the QCD evolution 
of parton densities, electroweak couplings and the prop- 
agator mass of the Z° gauge boson. Furthermore, such 
measurements allow the searches for physics beyond the 
Standard Model, such as resonance searches or contact 
interactions. 

2 Phenomenology of F2 at low x and low Q 2 

2. 1 Phenomenology of the low Q 2 region 

The primary purpose is to use NLO DGLAP QCD on the 
one hand and the simplest non-perturbative models on 
the other to explore the Q 2 transition region and through 
probing their limitations to shed light on how the pQCD 
description of F2 breaks down. One way to understand 
the rise in F 2 at low x is advocated by Gliick, Reya and 
Vogt (GRV94) who argue that the starting scale for the 
evolution of the parton densities should be very low (~ 
0.3 GeV 2 ) and at the starting scale the parton density 
functions should be non-singular. The observed rise in 
F 2 , with a parameterisation valid above Q 2 « 1 GeV 2 , is 
then generated dynamically. On the other hand, at low 
x one might expect that the standard NLO Q 2 evolution 
given by the DGLAP equations breaks down because of 
the large ln(l/x) terms that are not included. Such terms 
are taken into account by the BFKL formalism, which in 



leading order predicts a rising F2 at low x. The rise 
comes from a singular gluon density, xg ~ x A , with A 
in the range —0.3 to -0.5. Clearly accurate experimental 
results on F% at low x and the implied value of A are of 
great interest. 

At some low value of Q 2 pQCD will break down and 
non-perturbative models must be used to describe the 
data. At low x and large j*p centre-of-mass energy, W 
yj Q 2 /x, the total j*p cross-section is given by 

47T 2 a 



cr T + a L 



F 2 (x,Q 2 ) (1) 



where o~t and ol are the cross-sections for transversely 
and longitudinally polarised virtual photons respectively. 
Two non-perturbative approaches are considered, the 
generalised vector meson dominance model (GVMD) 
and a Regge-type two component Pomeron+Reggeon ap- 
proach a la Donnachie and Landshoff (DL) to give a good 
description of hadron-hadron and photoproduction total 
cross-section data. 

2.2 Measurement of F 2 with Shifted Vertex Data 

The shifted vertex data correspond to an integrated lu- 
minosity of 236 nb _1 taken in a special running period, 
in which the nominal interaction point was offset in the 
proton beam direction by +70 cm, away from the detect- 
ing calorimeter. Compared to the earlier shifted vertex 
analysis, for the 1995 data taking period the calorimeter 
modules above and below the beam were moved closer 
to the beam, thus extending the shifted vertex Q 2 range 
down to 0.6 GeV 2 . 

The double differential cross-section for single 
virtual-boson exchange in DIS is given by 



dx dQ 2 



a the negative of the square of the four-momentum transfer between 
the positron and the proton 
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Figure 1: Low Q 2 F2 data for different Q 2 bins together with the 
ZEUSDL style Regge model fit to the ZEUS BPC95 data. At larger 
Q 2 values the ZEUS NLO QCD fit is also shown. 



where R is related to the longitudinal structure function 
F L by R = F L /(F 2 - F L ) and S r gives the radiative 
corrections to the Born cross-section, which in this kine- 
matic region is at most 10%. The parity violating term 
xF% arising from the Z° exchange is negligible in the Q 2 
range of this analysis. Further details about the data 
analysis can be found in ref. u. 

Fig. 1 shows the results for F2 as a function of x 
in bins of Q 2 between 0.65 and 6 GeV 2 (ZEUS SVX95) 
together with ZEUS F2 measurements at very low Q 2 = 
0.11 - 0.65 GeV 2 (ZEUS BPC95) and at larger Q 2 those 
from the ZEUS94. There is good agreement between the 
different ZEUS data sets in the region of overlap. Also 
shown are data from the shifted vertex measurements by 
HI (HI SVX95) and fixed target data from E665. The 
steep increase of F2 at low x observed in the higher Q 2 
bins softens at the lower Q 2 values of this analysis. The 
curves shown will be discussed later in the text. 



2.3 The low Q 2 region 



We first give an overview of the low Q 2 region, Q 2 < 
5 GeV 2 , taking ZEUS SVX95, BPC95 and ZEUS94 F 2 
data. Using Eq. [l] we calculate a] Q f values from the 
F2 data. The DL model predicts that the cross-section 
rises slowly with energy oc W 2X , A = ap - 1 w 0.08 
and this behaviour seems to be followed by the data at 
very low Q 2 . Above Q 2 = 0.65 GeV 2 , the DL model 
predicts a shallower rise of the cross-section than the data 
exhibit. For Q 2 values of around 1 GeV 2 and above, the 



GRV94 curves describe the qualitative behaviour of the 
data, namely the increasing rise of aj a f with W 2 , as 
Q 2 increases. This suggests that the perturbative QCD 
calculations can account for a significant fraction of the 
cross-section at the larger Q 2 values. 

For the remainder of this section we concentrate 
on non-pcrturbative descriptions of the ZEUS BPC95 
data. Since BPC95 data are binned in Q 2 and y wc first 
rewrite the double differential cross-section of Eq. || as 
-JLj^ = r • (cr T + ea L ) where a L = jj^F L and a T has 
been defined by Eq. [l| The virtual photon has flux factor 
r and polarisation e. Keeping only the continuum states 
in the GVMD at a fixed W the longitudinal and trans- 
verse j*p cross-section are related to the corresponding 
photoproduction cross-section Oq P by 
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where the parameter £ is the ratio o~ v L v /o~ T 



^ p /cr^ p for vector 
meson (V) proton scattering and Mq is the effective vec- 
tor meson mass. Neither £ nor Mq are given by the 
model and they are either determined from a fit to data 
or by other approaches. As we do not have much sen- 
sitivity to £ and it is small (0.2 - 0.4) we set it here to 
zero. We thus have 9 parameters to be determined by 
fitting the BPC data to the simplified GVMD expression 

n 2 , ,,2 IP 

F2 = IA2 tV in 8 bins of W between 104 and 251 
GeV. The fit is reasonable and its quality might also be 
judged from the upper plot in Fig. 2. The value obtained 
for Af 2 is 0.53±0.04(stat)±0.09(sys). The resulting ex- 
trapolated values of <rj p are shown as a function of W 2 in 
the lower plot of Fig. 2, along with measurements from 
HERA and lower energy experiments. The extrapolated 
BPC data lie somewhat above the direct measurements 
from HERA. They are also above the cross section pre- 
diction of the DL model. It should be clearly understood 
that the ctJ p data derived from the BPC95 data are not a 
measurement of the total photoproduction cross-section 
but the result of a physically motivated ansatz. 

The simple GVMD approach just described gives a 
concise account of the Q 2 dependence of the BPC data 
but it says nothing about the energy dependence of <7q P . 
To explore this aspect of the data we turn to a two com- 
ponent Regge model 
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where P and R denote the Pomeron and Reggeon contri- 
butions. The Reggeon intercept an is fixed to the value 
0.5 which is compatible with the original DL value and by 
the re-evaluation of Cudell et al. With such an intercept 
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Figure 2: Upper plot: ZEUS BPC95 measurements of the total 
cross-section <tt + eox in bins of W and the GVMD fit to the data. 
Lower plot: a^ t as a function of W 2 . The ZEUS BPC95 points 
are those from the GVMD extrapolation of <Tq P . 



the Reggeon contribution is negligible at HERA energies. 
Fitting the extrapolated BPC95 data alone yields a value 
1.141 ±0.020(siai) for ap. Fitting both terms to the real 
photoproduction data (with W 2 > 3 GeV 2 ) and BPC95 
data yields ap — 1.101 ± 0.002(siai). Including in addi- 
tion the two original measurements from HERA as well 
gives ap = 1.100 ± 0.002(stai). All these values of ap 
are larger than the value of 1.08 used originally by DL, 
but we note that the best estimate of Cudell et al. is 
1.0964^0 0094) which within the errors is consistent with 
our result. The final step in the analysis of the BPC data 
is to combine the GVMD fitted Q 2 dependence with the 
Regge model energy dependence 
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The parameters Mq and an are fixed to their previous 
values of 0.53 and 0.5, respectively. The 3 remaining pa- 
rameters are determined by fitting to real photoproduc- 
tion data and the original BPC data. The description of 
the low Q 2 F 2 data given by this DL style model is shown 
in Fig. 1. Data in the BPC region Q 2 < 0.65 GeV 2 is 
well described. At larger Q 2 values the curves fall below 
the data. Also shown in Fig. 1 for Q 2 > 6 GeV 2 are 
the results of a NLO QCD fit (full line) as described in 
Sec. |J. 

24 F 2 slopes: d\nF 2 /d\n(l/x);dF 2 /d\nQ 2 

To quantify the behaviour of F 2 as a function of 
Q 2 and x at low x we calculate the two slopes 



dlnF 2 /d]n(l/x);dF2/dlnQ 2 from the ZEUS SVX95, 
BPC95 and ZEUS94 data sets. 

At a fixed value of Q 2 and at small x the behaviour 
of F 2 can be characterised by F 2 oc x~ A , with A tak- 
ing rather different values in the Regge and BFKL ap- 
proaches. X e ff is calculated from horizontal slices of 
ZEUS F 2 data between the y = 1 HERA kinematic limit 
and a fixed cut of x < 0.01, here including E665 data. 
In a given Q 2 bin ( x ) is calculated from the mean value 
of ln(l/a;) weighted by the statistical errors of the cor- 
responding F 2 values. The same procedure is applied to 
the theoretical curves shown for comparison. 

Figure 3 shows the measured values of X e f / as a func- 
tion of Q 2 . From the Regge approach one would ex- 
pect X e ff « 0.1 and independent of Q 2 . Data for Q 2 < 
1 GeV 2 is consistent with this expectation. The linked 
points labelled DL are calculated from the Donnachie- 
Landshoff fit and as expected from the discussion of 
the previous section are somewhat below the data. For 
Q 2 > 1 GeV 2 , X e ff increases slowly to around 0.3 at 
Q 2 values of 40 GeV 2 . Qualitatively the tendency of 
Xeff to increase with Q 2 is described by a number of 
pQCD approaches. The linked points labelled GRV94 
are calculated from the NLO QCD GRV94 fit. Although 
the GRV94 prediction follows the trend of the data it 
tends to lie above the data, particularly in the Q 2 range 
3 -20 GeV 2 . 
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Figure 3: Left plot: din F2 / dln(l/ x) as a function of Q 2 calculated 
by fitting ZEUS and E665 F 2 data in bins of Q 2 . Right plot: 
dF2 /d In Q 2 as a function of x calculated by fitting ZEUS F2 data 
in bins of x. 

Within the framework of pQCD, at small x the be- 
haviour of F 2 is largely determined by the behaviour of 
the sea quarks F 2 ~ xS, whereas the dF 2 /dlnQ 2 is de- 
termined by the convolution of the splitting function P qg 
and the gluon density, dF 2 /dlnQ 2 oc a s P qg ®g. In order 
to study the scaling violations of F 2 in more detail the 
logarithmic slope dF 2 jd In Q 2 is derived from the data by 
fitting F 2 — a + 61nQ 2 in bins of fixed x. The statistical 
and systematic errors are determined as described above. 
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The results for (IF2/ 'd\nQ 2 as a function of x are shown 
in Fig. 3. For values of x down to 3 x 10 -4 , the slopes 
are increasing as x decreases. At lower values of x and 
Q 2 , the slope decreases. Comparing the rapid increase 
in F2 at small x with the behaviour of the dF^j 'd\nQ 2 , 
one is tempted to the naive conclusion that the under- 
lying behaviour of the sea quark and gluon momentum 
distributions must be different at small x, with the sea 
dominant and the gluon tending to zero. The failure of 
DL is in line with the earlier discussion. GRV94 does not 
follow the trend of the data when it turns over. 

2.5 NLO QCD fit to F 2 data 

In perturbative QCD the scaling violations of the F2 
structure function are caused by gluon bremsstrahlung 
from quarks and quark pair creation from gluons. In 
the low x domain accessible at HERA the latter process 
dominates the scaling violations. A QCD analysis of F2 
structure functions measured at HERA therefore allows 
one to extract the gluon momentum density in the proton 
down to low values of x. In this section we present NLO 
QCD fits to the ZEUS 1994 nominal vertex data and the 
SVX95 data of this paper. We are not attempting to 
include all available information on parton densities, but 
concentrating on what ZEUS data and their errors allow 
us to conclude about the gluon density at low x. 

To constrain the fits at high x proton and deuteron 
F 2 structure function data from NMC and BCDMS are 
included. The kinematic range covered in this analysis is 
3 x KT 5 < x < 0.7 and 1 < Q 2 < 5000 GeV 2 . 

The QCD predictions for the F2 structure functions 
are obtained by solving the DGLAP evolution equations 
at NLO. These equations yield the quark and gluon mo- 
mentum distributions at all values of Q 2 provided they 
are given at some input scale Q 2 . In this analysis we 
adopt the so-called fixed flavour number scheme where 
only three light flavours (u, d, s) contribute to the quark 
density in the proton. The corresponding structure func- 
tions F£ and F| are calculated from the photon-gluon 
fusion process including massive NLO corrections. The 
input valence distributions are taken from the parton 
distribution set MRS(R2). As for MRS(R2) we assume 
that the strange quark distribution is a given fraction 
K s = 0.2 of the sea at the scale Q 2 = 1 GeV 2 . The 
gluon normalisation is fixed by the momentum sum rule. 
The input value for the strong coupling constant is set 
to a s (M|) = 0.118 and the charm mass is taken to be 
m c = 1.5 GeV. In the QCD evolutions and the evalua- 
tion of the structure functions the renormalisation scale 
and mass factorisation scale are both set equal to Q 2 . In 
the definition of the \ 2 only statistical error are included 
and the relative normalisation of the data sets is fixed 
at unity. The fit yields a good description of the data as 



shown in Fig. 1. We have also checked that the gluon ob- 
tained from this fit to scaling violations is in agreement 
with the recent ZEUS measurements of charm produc- 
tion and F£ in deep inelastic scattering at HERA. 

Two types of systematic uncertainties have been con- 
sidered in this analysis. 'HERA standard errors' contain 
statistical error on the data, experimental systematic un- 
certainties, relative normalisation of the different data 
sets and uncertainties on a s , the strange quark content 
of the proton and the charm mass. 'Parametrisation er- 
rors' contain uncertainties from a x 2 definition including 
statistical and experimental systematic errors, variations 
of the starting scale Q 2 , and an alternative, more flexible 
parametrisation of the gluon density using Chebycheff 
polynomials. The first type of errors amounts to 16% 
Ag/g at x = 5 x 10~ 5 , Q 2 = 7 GeV 2 , the second type 
yields 9.5% in Ag/g. 
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Figure 4: The quark singlet momentum distribution, xT, (shaded), 
and the gluon momentum distribution, xg(x) (hatched), as a func- 
tion of x at fixed values of Q 2 = 1, 7 and 20 GeV 2 . The error bands 
correspond to the quadratic sum of all error sources considered for 
each parton density. 



The three plots of Fig. 4 show the distribution for 
and xg as a function of x for Q 2 at 1, 7 and 20 GeV 2 . 
It can be seen that even at the smallest Q 2 xT, is rising 
at small x whereas the gluon distribution has become 
almost flat. These results give support to the naive con- 
clusion of Sec. 2.4, that the sea distribution dominates 
at low x and Q 2 . At Q 2 = 1 GeV 2 the gluon distribution 



is poorly determined and can, within errors, be negative 
at low x. 
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3 Measurement of the Proton Structure Func- 
tion F 2 from 1996 and 1997 data 

3.1 Kinematics in Deep Inelastic Scattering 

Recalling the double differential NC cross-section (||) , but 
now including the corrections (5l and S3) for and xF 3 
yields 



d 2 a 
dx dQ 2 



2na 2 Y+ 



F 2 (1 - 6 L - 6 3 ) (1 + 6 r ) (5) 



Here the F 2 structure function contains contributions 
from virtual photon and Z° exchange 
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where Mz is the mass of the Z° and Fj 
are the contributions to F 2 due to photon exchange, Z° 
exchange and jZ° interference respectively. In this anal- 
ysis we determined the structure function F 2 m using 1996 
and 1997 data with an integrated luminosity of 6.8 pb _1 
and 27.4 pb _1 , respectively. 

The selection and kinematic reconstruction of NC 
DIS events is based on an observed positron J,nd the 
hadronic final state. For further details see ref. □. 

3.2 Results 

Monte Carlo samples are used to estimate the accep- 
tance, migration, radiative corrections, electroweak cor- 
rections and background contributions. F 2 m is then de- 
termined based on a bin-by-bin unfolding. The result- 
ing statistical error, including the Monte Carlo statis- 
tics, ranges from 2% below Q 2 = 100 GeV 2 to 5-6% at 
Q 2 «800 GeV 2 . 

The systematic uncertainties have been estimated by 
varying the selection cuts, efficiencies and reconstruction 
techniques and redetermining the cross section includ- 
ing background estimates. Potential error source such as 
possible detector misalignment, event vertex reconstruc- 
tion, calorimeter energy scale, positron identification ef- 
ficiency, background contributions and hadronic energy 
flow have been considered. The total systematic uncer- 
tainty amounts to 3-4% except at low and high y, where 
it grows to 12%. At the present preliminary state of the 
analysis we estimate an overall normalisation uncertainty 
of 3%. 

The resulting F 2 m is shown as a function of x for 
fixed Q 2 in Figure 5. Results from our previous anal- 
ysis, and from fixed target experiments are also shown 
for comparison. At low Q 2 the rise in F 2 for x — > is 
measured with improved precision. The coverage in x 
has also been extended to higher x, yielding extended 



overlap with the fixed target experiments; in the over- 
lap region reasonable agreement has been found. The F 2 
scaling violation from this analysis and the fixed target 
data are also shown in Figure 5. For Q 2 > 100 GeV 2 
the increase in statistics allows a measurement of F 2 m in 
smaller bins with respect to our previous measurement. 
Above Q 2 = 800 GeV 2 , the statistical error grows typ- 
ically to 5-15% and dominates the total error. Overall 
our data are in agreement with our published data. 
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Figure 5: Top and bottom left plots: F^ m versus x for fixed Q 
Bottom right plot: FS m as a function of Q 2 for fixed x. 
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